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Arrays of atoms trapped in optical lattices are appealing as storage media for photons, since 
motional dephasing of the atoms is eliminated. The regular lattice is also associated with band 
structure in the dispersion experienced by incident photons. Here we study the influence of this band 
structure on the efficiency of quantum memories based on electromagnetically induced transparency 
(EIT) and on Raman absorption. We observe a number of interesting effects, such as both reduced 
and superluminal group velocities, enhanced atom-photon coupling and anomalous transmission. 
These effects are ultimately deleterious to the memory efficiency, but they are easily avoided by 
tuning the optical fields away from the band edges. 

PACS numbers: 42.50.Ex, 42.50.Ct, 42.50.-p 
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I. INTRODUCTION 

Quantum memories for photons based on coherent ab- 
sorption in atomic ensembles promise to provide the 
Ught-matter interface required for scalable quantum op- 
tical networking [IHS] . Recently, arrays of atoms trapped 
in optical lattices [3] have been used as a storage medium 
[II [8] . The lattice eliminates decoherence via atomic dif- 
fusion and collisions, and coherence times of many sec- 
onds are feasible [S]. On the other hand, the periodic 
arrangement of atoms is expected to generate a photonic 
band structure for incident signals, with certain optical 
frequencies forbidden due to the destructive interference 
of scattered fields [TUHH]- At frequencies close to the 
edge of a forbidden band, the group velocity of an inci- 
dent signal is reduced, and the signal therefore interacts 
with the atoms for longer [Tni IS] . It is not clear what 
this implies for the efficiency of optical lattice quantum 
memories compared with their free-atom counterparts. 
In this paper we model quantum storage in a periodically 
structured ensemble in order to investigate this question. 
Our results indicate that while the interaction of the sig- 
nal with the atoms grows stronger near a band edge, the 
memory efficiency is reduced, because the coupling to the 
storage state of the memory is not enhanced. While the 
lattice band structure is not advantageous for quantum 
memory, simply tuning the optical fields away from the 
band edges allows for efficient storage, with all the ben- 
efits of reduced decoherence that accrue from the ability 
to fix the atoms in space. 



II. MODEL 

Optical lattices are capable of supporting highly entan- 
gled quantum states, ranging from atomic superfluids to 
Mott insulators [6'. However, a semiclassical treatment 
in which the optical fields and the atomic positions re- 
main un-quantized is sufficient to analyze the efficiency of 
quantum storage, which quantity involves only normally 
ordered products of field amplitudes. The following anal- 
ysis therefore does not depend on the spatial correlations 



within the lattice, and so it also applies to any ensemble 
memory in which the atomic density varies periodically 
in space on the scale of an optical wavelength, such as 
might be produced by doping a photonic crystal structure 

We consider quantum memories based on EIT j31 ITl] 
and Raman scattering |3l|T5l[T6] in an ensemble of A- type 
atoms (see figure fl]). In these memory protocols, the sig- 
nal photon to be stored is absorbed on the 1 1) -^ | 2) 
transition, and mapped into the storage state | 3) by an 
intense control field. In a Raman memory the signal and 
control are tuned far from resonance by a common detun- 
ing A; the fields are tuned into resonance, with A = 0, 
in an EIT memory. A single theoretical model therefore 
suffices to describe both protocols [17l [18] . We consider 




FIG. 1: (Color online) Quantum memory in an optical lat- 
tice, (a) A signal field A is directed into an ensemble of atoms 
trapped in a regular array, along with a bright control field 
n. (b) Each atom has a A-Ievel structure. The control and 
signal fields are tuned into two-photon resonance, precipitat- 
ing the transfer of atoms from the the ground state | 1) to 
the long-lived storage state | 3) . In an EIT memory, the com- 
mon detuning A of the signal and control from the excited 
state I 2) is zero; in a Raman memory, A S> 7, where 7 is the 
homogeneous linewidth of the | 1) -f-)- | 2) transition. 



propagation of the signal in one dimension, along the z- 
axis. The signal electric field Eg couples to the atomic 
polarization P, via the wave equation 
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E, = -y.Qd'lP, 



(1) 



Bragg scattering of the signal from the periodic array of 
atoms in the optical lattice introduces backward travel- 
ling components into the signal beam. We account for 
this by introducing a carrier wave for the signal 0, which 
is not in general a plane wave. We introduce the slowly 
varying amplitude A for the signal by the relation 



Es{z,t) = igsA{z,t)(p{z)e 



(2) 



where gs = ^JhujgjIei^Ac^ with A the cross-sectional area 
of the signal field, and where uJs is the signal carrier 
frequency. We describe the modulation of the atomic 
density due to the lattice with the periodic function 
m(z) — m{z -\- a), where a is the lattice constant, and 
where /„ m{z) dz — L, with L the length of the ensem- 
ble, so that m — 1 describes a uniform ensemble with no 
optical lattice. The slower spatial variation of the atomic 
polarization over the length of the ensemble is described 
by the amplitude P, which we define by the relation 
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-m,{z)P{z, t)(f){z)e 



(3) 



where we have factorized out both the lattice modulation, 
and the signal carrier wave (j). Here n is the constant 
average number density of atoms in the ensemble, and di2 
is the dipolc matrix element for the | 1) o | 2) transition. 
The atomic dynamics in the presence of the signal and 
control fields are described by the Bloch equations, which 
can be written in the form [T7l [18] 



dtP 
dtB 
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ifl*P. 



ikA + iflB, 



(4) 
(5) 



where F = 7 — iA is the complex detuning, with 7 
the homogeneous linewidth of the excited state 1 2) , 
and where the coupling constant is k = y/drf/L, with 
d — di2^snL /2eQch'y the resonant optical depth of the 
ensemble. B represents the amplitude of the long-lived 
1 1) O I 3) Raman coherence, or spin wave, into which the 
signal field is mapped by the memory interaction. The 
slowly varying Rabi frequency fl ~ fl(t— z/c) represents 
the temporal profile of the control pulse, which experi- 
ences no dispersion and travels at c, since it couples the 
states I 2) and | 3) , whose populations remain negligible 
at all times. 

When the spectral bandwidth 5 of the control field is 
sufficiently narrow, S <^ A, or S <^ dj, the atomic polar- 
ization P can be adiabatically eliminated |191l20j . This is 
achieved by setting the time derivative on the left hand 
side of Eq. Q to zero, and substituting the resulting solu- 
tion for P into Eqs. ([5) and ([l]), using Eq. ^. We make 
a slowly varying envelope approximation by neglecting 
the terms d'^A, d^A and both dtP and d^P in Eq. ^. 
This yields the pair of coupled equations 

2d,Ad,(l) + Ad^,<j> 

2i/c, 
(l>dtA 
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2ksmn — B(j), 



B = iVL*-A, 



(6) 

(7) 



where V — V{z) = 1 + ^PL™ ' with kg = uis/c. These 
equations simplify considerably if we require that the car- 
rier wave (j) should satisfy the Schrodinger-like equation 



[dl + klV{z)] <j>{z) = 0. 



(8) 



Here V plays the role of a potential, which takes the form 
of a constant with a small modulation added to it. The 
periodicity of the potential ensures that the carrier wave 
can be expressed in the Bloch-Floquet form [21] 
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where k is the crystal momentum and where Ukv{z) = 
Uku{z + a) is a periodic Bloch function. Substituting 
Eq. ([9]) into Eq. ((Sl), we find that the band index 1/ must 
be chosen along with the crystal momentum k in order 
to satisfy the equation 

[af + 2ikd, - fc2 + klV{z)] UkAz) = 0. (10) 

If it is not possible to find values of k and v that satisfy 
Eq. (10), the signal frequency Wg is said to lie within a 



photonic band gap. We will be interested in frequencies 
that lie close to the edge of a band gap. 

When Eq. ( 10 1 is satisfied, Eq. ^ becomes 

— (j}dtA = ksmn-B(t). (11) 
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The potential V contains the complex detuning, whose 
imaginary part describes absorption of the signal field by 
the atoms. In general, therefore, the crystal momentum 
k has a small but non-zero imaginary part, and the car- 
rier wave is damped to some degree. It is convenient, in 
analyzing the memory efficiency, to transfer this damp- 
ing to the slowly varying signal amplitude A, which we 
do by making the transformations A — > Ae~^'^^''^^ and 

A final simplification is achieved by 'projecting' 
Eq. (11) onto the mode (j) [13]. This is accomplished by 



multiplying both sides of Eq. (11) by a conjugate mode 
-0, and integrating over a unit cell of the lattice, using 
the fact that the functions A, B and ft are effectively 
constant over this range. The conjugate mode ijj is not 
equal to (j)* , because Eq. (l8| is not generally a Hermitian 
eigenvalue equation. Re- writing Eq. ([8]), we have that 
is the right eigenvector of the operator M = —V~^&1, 
with eigenvalue fc^. The corresponding left eigenvector 
of M is then equal to ip- 

After these manipulations, we arrive at the following 
set of equations of motion for the quantum memory 

[a,+Im{fc}]A = -{c/v^)a-^n{T + Pz)%2) 



dr 



|r!(T + /3z)|^ 



B 



?0*(r + /?z)A, 



(13) 



where t = t — z/v^ denotes the time in a frame moving 
at the signal group velocity, given by [13] 



i^{z) [~\dz(j){z)\ dz, 



(14) 



and where /? = 1/ug — 1/c quantifies the rate at which 
the control walks off from the signal [22] . The overlap a 
captures the degree to which the optical Bloch mode is 
matched to the lattice modulation 



ijj(z)(f)(z)m(z) dz 



(15) 



Note that far from a band edge, the carrier wave (p takes 
the form of a plane wave, so that a = 1, Vg = c and (3 ~ 
0. The equations then reduce to the standard equations 
describing the Raman and EIT protocols in a disordered 
ensemble P^ITT]. 



III. EFFICIENCY 

To analyze the memory efficiency, we need to solve the 
coupled system Eqs. ([l2]), ( [l3| and construct the Green's 
function K such that^22j 



Boutiz) 



K{z,T)AiniT)dT, 



(16) 



where Bout{z) — B(z, r — >■ oo) is the spin wave left in the 
atoms at the end of the storage interaction, and where 
An(''') = A{z = 0,t) is the amplitude of the signal im- 
pinging on the entrance face of the ensemble. The stor- 
age efficiency is defined as the ratio of the number of final 
spin wave excitations to the number of incident photons, 
which is given by 



^ /J-|i3o„t(z)pdz 
^ JZlAUrWdr- 



(17) 



The optimal storage efficiency rjopt can be found from the 
singular value decomposition [23, of the Green's function 
K by squaring the largest singular value [24]. 

In general, the equations of motion cannot be solved 
analytically, because of the walk-off between the signal 
and the control field. However, the typical size of optical 
lattices (~ 1 mm) is small compared with the longitudi- 
nal spatial extent Tc of a typical photonic wavepackct, 
where T is the signal pulse duration, for all but the short- 
est signal pulses. In this case, even very close to a band 
edge, we have f3L <^ T, and we can safely drop /3 from 
Eqs. ( [l2| ), (13 1. The Green's function is then given by 



K{z,t) = ^^(T)e-'^(^'^) Jo 2JaKzuj{T)c/vg/r 



(18) 
where xiz,T) — Im{/i;}z 4^ lu{t)/T, and where we 
have defined the integrated Rabi frequency a;(r) — 
J^^ |^('''')l dr'. Here Jo denotes the zero'th order ordi- 
nary Bessel function of the first kind. A coordinate trans- 
formation from T to w = w(t) reveals that the singular 
values of this kernel do not depend on the temporal pro- 
file J7; they depend only upon the energy in the control 



pulse, represented by the limit u!{t —>■ oo). Therefore the 
optimal storage efficiency for a lattice memory depends 
on the control pulse energy, but not on its shape. 

To understand the effect of the photonic band struc- 
ture induced by the optical lattice, we study the varia- 
tion of the optimal efficiency yyopt as we vary the lattice 
constant a, with all other quantities fixed. In the cal- 
culations, we assume atomic parameters similar to those 
for cesium or rubidium, so we set As — 271 /ks — 800 nm, 
and T2 — 1/7 = 30 ns. We assume that the atoms are 
trapped in an optical lattice of length L — \ mm, and of 
width sufficient to cover the full beam waist of the signal 
field. We assume a Gaussian profile m{z) oc e"*^^/*^'^ for 
the density modulation due to the lattice potential, with 
width w = a/10. For concreteness, a Gaussian profile 
^{t) = r^oe^'^'^/^^ is assumed for the control pulse, with 
a peak Rabi frequency JIq = 5.5/T. 

We consider two situations. First, a broadband off- 
resonant Raman protocol, with d = 300, T = 3 ns and 
A = 15/T, and second, a narrowband resonant EIT pro- 
tocol with d = 30, T = 30 ns and A = [I8]. Note that 
the adiabatic condition Td-y ^ 1 is well-satisfied in both 
cases [T3- Both types of memory exhibit band structure 
(Figure ^ , but the dispersive properties of Raman and 
EIT memories are quite different. 

Figure [3] shows the variation of the optimal storage 
efficiency yyopt, along with the group velocity Vg and the 
overlap parameter a, for the two memory protocols, as 
the frequency of the signal field approaches a band edge. 
In the Raman memory, the group velocity falls steadily 
to zero, and the overlap increases, consistent with the 
transition via Bragg scattering of the signal carrier from 
a plane to a standing wave with nodes between the lattice 
sites. Both of these effects might be expected to increase 
the memory efficiency, but instead the memory efficiency 
actually falls. 

The reason for this surprising negative result appears 
to be related to the joint atomic/optical character of the 
memory interaction: although close to a band edge the 
coupling of the signal field to the atomic polarization in- 
creases, the coupling of the signal field to the storage 
state I 3) is unchanged, since this is mediated by the con- 
trol field, which experiences no dispersion. The atomic 
polarization simply re-radiates the signal field at an en- 
hanced rate, and less of the signal is mapped into the 
desired storage state. 

The behaviour of the EIT memory is quite different. 
Here, the modulation of the potential V{z) is purely 
imaginary, since A = 0. This gives the band structure a 
radically different appearance (see part (b) of Figure l2| , 
and this is manifested in the variation of the group ve- 
locity, which becomes superluminal briefiy, before falling 
sharply, as shown in part (d) of Figure pi The overlap 
parameter also falls, in contrast to the Raman case, sug- 
gesting that at the band edge the EIT carrier becomes a 
standing wave with nodes aligned with the lattice sites. 
The optimal efficiency changes rather counter-intuitively. 
It first decreases, and then rises sharply. We understand 
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FIG. 2: (Color online) Dispersion in an optical lattice, (a) 
The dispersion relation for the signal field: that is, the de- 
pendence of the Bloch wave eigenvalue ks upon the real part 
of the crystal momentum k. The dispersion relation is multi- 
valued; for a given crystal momentum, the eigenvalues are 
enumerated by the band index i/. The bands with v — 1 
and J/ = 2 are shown here. The band gaps appear at the 
position indicated by the arrow, which marks the edge of the 
first Brillouin zone, (b) Close-up of the first Brillouin zone 
edge for the case of a Raman memory. The band gap has 
the appearance of a 'standard' anti-crossing, (c) Correspond- 
ing close-up for an EIT memory. The band structure here is 
more complex, because the modulation on the potential V{z) 
is imaginary. The gradient of the v — 1 band increases at 
first, before sharply approaching zero. The bandgap itself is 
too small to see in the plot. 



this behaviour in the following way. 

In an EIT memory, the coupling of the signal to the 
atomic polarization is directly related to the memory effi- 
ciency; the control field dresses the atoms and mixes the 
storage and excited states, excitation of the atomic po- 
larization then directly excites the storage state. There- 
fore the initial reduction in a and the increase in Vg 
act to decrease the memory efficiency. The subsequent 
'turnaround' in efficiency very near to the band edge 
is partly explained by the sharp reduction in Vg, which 
serves to increase the atom-signal coupling. But a second 
effect is more significant: anomalous transmission. This 
refers to a reduction in the absorption that is character- 
istic of an imaginary potential, when close to a forbidden 
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FIG. 3: (Color online) Storage efficiency, group velocity and 
overlap near a band gap. (a) The storage efficiency of the 
Raman memory is plotted as a function of the signal frequency 
near the band gap. The blue solid line is found from the 
analytic kernel in Eq. (18 1; the black crosses are found from 
a numerical solution of Eqs. (121, ( |f3[ ), which takes account 
of walk-off between the signal and control. The efficiency 
falls as the signal frequency approaches the band edge, (b) 
The corresponding variation of both the group velocity, which 
falls, and the overlap parameter a, which increases as the 
band edge is approached, (c) EIT storage efficiency, which 
falls, and then rises sharply as the signal field approaches the 
band edge. In part (d) we plot the group velocity, which rises 
and then falls, and the overlap parameter a, which falls, as 
the band edge is approached. 



band — it has been observed in the transmission of atoms 
through an optical standing wave p5H?7] ; here it pertains 
to the transmission of an optical field through a periodic 
array of atoms. The attenuation of the signal carrier for 
the Raman and EIT protocols is shown in Figure |4j In 
the absence of any modulation, m{z) = 1, and Eq. (|8]) 



can be solved trivially, yielding Ini{fc} = Re{d7/rL}. 
Any departure from this value is associated with the 
band structure arising from the lattice modulation. We 
quantify these departures by defining the dimensionless 
damping parameter /i = 1 — lia {k} L /He {d'j /T} . The 
variation of fi as the signal frequency approaches a band 
edge is plotted in Figure |4] for both memory protocols. 
In the Raman case, the damping increases near the for- 
bidden band, consistent with the transition of the signal 
carrier from a propagating to an evanescent wave. In 
the case of EIT however, the damping drops, becoming 
negative near the band gap. The absorptive scattering 
of the signal is thus dramatically reduced, and it is this 
anomalous transmission [25H27] that explains the sudden 
increase in the optimal storage efficiency when very near 
the band edge. 

Here we note that our theoretical model has limited 
applicability so close to the band edge. The slowly vary- 
ing envelope approximation used in deriving Eqs. (rol, ([7]) 
requires that a single carrier wave correctly describes the 
propagation of all the frequencies comprising the signal 
pulse. When the dynamics become extremely dispersive, 
as they do in very close proximity to the stop band, the 
bandwidth of the signal pulse will itself span a range of 
carrier modes. The detailed analysis of this effect lies 
beyond the scope of the present treatment, but in any 
case it will lead to the dispersive break-up of the signal 
pulse, and to a reduced memory efficiency. Therefore one 
should probably not conclude that enhanced storage ef- 
ficiency can be achieved via EIT by tuning very close to 
a band edge. In the next section we show that reflection 
losses negate any such enhancement in any case. 




FIG. 4: (Color online) Damping, (a) The variation of the 
damping parameter fi = 1 — Im {k} L/Re {dy/T} for the Ra- 
man memory, as a function of the signal frequency, approach- 
ing the bandgap. In the absence of band structure, ^ — 0, 
but the increase in fi close to the band edge is associated with 
increased absorption of the signal carrier wave, (b) The cor- 
responding variation of n for the EIT memory. Close to the 
band edge, fi becomes increasingly negative, which is associ- 
ated with reduced absorption. This phenomenon is known as 
anomalous transmission [25H27| . 



IV. REFLECTION 

A final important consideration for the lattice memory 
is the refiectivity of the interface between free space and 



the ensemble [Tni[2Hl[2n]- Although the optimal efficiency 
of EIT appears to approach unity very near the band 
edge, it seems that the portion of the signal field lost 
due to refiection at the entrance face of the memory is 
sufficiently large to abrogate any advantage. 

We obtain an expression for the reflectivity of the en- 
semble by imposing continuity of the signal carrier waves 
and their derivatives at the free-space/lattice interface 
(see part (a) of Figure Isl) . The result is 
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n -r2 



n +r2 



where n = ksUki,{0), and r2 = kuk,^{0) - idzUku{z)\z=o- 
Parts (b) and (c) of Figure [s] show the dependence of the 
reflectivity, along with the appropriately attenuated stor- 
age efficiency (1 — R)riopt, on the proximity of the signal 
frequency to the band edge for the two memories. The re- 




3.12 3.125 3.13 3.135 3.14 2.8 2.85 2.9 2.95 3 3.05 3.1 
KoOj KcQj 



FIG. 5: (Color online) Reflection from the entrance face, (a) 
continuity of the incident, reflected and transmitted signal 
carrier waves at the entrance face of the lattice determine the 
reflectivity R = |rp of the memory, (b) the variation of R, 
and the corresponding storage efliciency (f — R)r]opt, as the 
signal frequency approaches the band edge, for the Raman 
memory, (c) the reflectivity and resulting efficiency of the EIT 
protocol. Both memory protocols suffer a dramatic reduction 
in efficiency close to the band edge, due to a sharp rise in the 
reflectivity of the lattice. 



fiectivity R increases dramatically as the signal frequency 
approaches the band gap, for both the EIT and Raman 
memory protocols. The fraction of the signal that pene- 
trates the ensemble and is subsequently stored is there- 
fore reduced, and the overall memory efficiency suffers 
at the band edge. In the case of the Raman memory, 
this effect compounds the reduction in efficiency arising 
from the dispersive propagation. For the EIT memory, 
the sharp rise in efficiency due to anomalous transmis- 
sion is tempered by the reflection losses, so that there is 



no longer any advantage in tuning the signal close to a 
band edge. 



V. CONCLUSION 

We have analyzed the effects of band structure on the 
efficiency of ensemble-based memory protocols in optical 
lattices. Although a number of interesting effects, includ- 
ing enhanced atom-light coupling, sub and superlumi- 
nal group velocities and anomalous transmission emerge 
from the model, the memory efficiencies are nonetheless 
reduced near a band edge. On the basis of this analy- 
sis, we conclude that one should avoid tuning the optical 
fields too close to a band edge when performing light 
storage in an optical lattice. Since a deep lattice po- 
tential is achieved by tuning the trapping lasers close to 
resonance, there is a real possibility of 'accidental' coin- 
cidence of the band edges with the signal frequency; this 



should be avoided if possible. 

Far from the forbidden bands, the memory efficiency 
in a lattice is the same as would be achieved in an equiv- 
alent disordered ensemble with the same optical depth. 
Of course, the long atomic coherence times available in 
optical lattices make them one of the most promising 
technological routes to useful quantum optical memories. 
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